Simulation of a spin-wave instability from atomistic spin dynamics 
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We study the spin dynamics of a Heisenberg model at finite temperature in the presence of an 
external field or a uniaxial anisotropy. For the case of the uniaxial anisotropy our simulations show 
that the macro moment picture breaks down. An effect which we refer to as a spin-wave insta- 
bility (SWI) results in a non-dissipative Bloch-Bloembergen type relaxation of the macro moment 
where the size of the macro moment changes, and can even be made to disappear. This relaxation 
mechanism is studied in detail by means of atomistic spin dynamics simulations. 

PACS numbers: 75.10.-b, 75.20.En, 75.40.Gb 



I. INTRODUCTION 

Relaxation processes for magnetization dynamics are 
poorly understood but play a crucial role for spin dy- 
namics in general. In this article we address one of 
the most fundamental processes of magnetization dy- 
namics, the uniform motion of the magnetization in an 
anisotropy field. The uniform motion of the magneti- 
zation and the relaxation of the uniform motion is of 
central importance in applications (magnetic switching 
in storage media etc.)ji We study here a ferromagnetic 
system which initially is excited by a finite angle rotation 
of the magnetization with respect to the anisotropy axis. 
This excitation brings the magnetization into a uniform 
motion where the magnetization eventually relaxes back 
to an alignment with the magnetization direction par- 
allel to the anisotropy axis. Different phenomenological 
models, such as Gilbert damping and Bloch-Bloembergen 
dampingf^ have been used for describing this macro 
level relaxation of the magnetization. In this article we 
perform simulations of magnetization dynamics on an 
atomic scale and we study the consequences for the macro 
scale behavior of the magnetization dynamics. 

The initial rotation of the magnetization of a ferro- 
magnet in an external field can be seen as an excitation 
of a large number of uniform k=0 magnons. During the 
relaxation process these magnons interact, dissipating en- 
ergy and angular momentum. Relaxation can occur via 
two processes, one where both energy and angular mo- 
mentum are transfered out (or in) of the magnetic sys- 
tem and the second where energy is transfered within the 
magnetic system, to other non-uniform k ^ magnons. 
The first process, which describes a dissipative damping 
in the equations of motion for magnetization dynamics, 
results in a Gilbert like relaxation, 
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where the effective field is assumed to lie in the z- 
direction and T is a relaxation parameter. This second 
process is the focus of this work. 

We address here a mechanism for the relaxation of 
the uniform motion of the magnetization within the 
spin-system itself which is seen to result in a Bloch- 
Bloembergen like damping of the magnetization of 
the system. Several such mechanisms exist, such as 
the Suhl instability, 4 2-magnon scattering^ 4-magnon 
scattering^ etc. All of the mentioned mechanisms rely 
on the dipolar interactions resulting in an energy lower- 
ing of non-uniform magnons and an energy degeneracy 
between the uniform magnons and certain non-uniform 
magnons. 

In this article we study a different mechanism, since 
dipolar interaction are not even included in our simula- 
tions. The mechanism that we address is instead due 
to the thermal fluctuations on an atomic scale of the 
magnetization combined with the nature of the uniax- 
ial anisotropy. Such a mechanism, which does not rely 
on dipolar interactions, were studied by Safonov et al& 
and recently by Kashuba£ and Garanin et a/i 10 ' 11 where 
it was shown that a spin wave instability (SWI) devel- 
ops in a uniaxial anisotropy field. As we show in this 
article, based on theoretical considerations, the insta- 
bility should develop on the atomic length scale as well 
as on the micrometer length scale which was treated by 
Kashuba, provided certain conditions are fulfilled. The 
instability is shown to be caused by the altering of the 
non-uniform thermal magnetic excitations of the system 
as it undergoes a uniform rotation. 
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where M is the macro moment, H the effective field, 7 
the gyromagnetic ratio, and a a damping parameter. The 
second process, which is described by the precessional 
term in the equations of motion, results in a special case 



II. DETAILS OF THE SPIN-DYNAMICS 
SIMULATIONS 

There are at least two approaches for studying mag- 
netization dynamics in simulations. Most common is 
what is pursued in micromagnetics, the solution of the 
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phenomenological Landau Lifshitz Gilbert (LLG) equa- 
tion on a micrometer length scale for a continuum 
magnetization.— An alternative approach, which is uti- 
lized here, is based on solving the equations of motion 
for magnetization dynamics where the magnetization on 
a nanoparticle or atomic scale is represented by a Heisen- 
berg Hamiltonian. The Heisenberg Hamiltonian is often 
successful in describing magnetic systems on an atomic 
scale, especially when using first principles calculations of 
the interatomic exchange. The current method is based 
on atomistic spin dynamics which has as starting point 
a quantum mechanical description from density func- 
tional theory of the evolution of the atomic spins. Other 
works which have taken this approach can be found in 
Refs. [IH1IS 

Our simulations are performed using the ASD (Atomic 
Spin Dynamics) package^ which is based on an atom- 
istic approach of spin dynamics. Interatomic exchange 
and magneto-crystalline anisotropy (MA) are included 
in the Hamiltonian. We use a parameterization of the 
interatomic exchange part of the form of a Heisenberg 
Hamiltonian, where the exchange parameters are calcu- 
lated from first principles theory. The effect of temper- 
ature is modeled by Langevin dynamics. Connection to 
an external thermal bath is modeled with a Gilbert like 
damping. The simulations are performed on bcc Fe using 
four coordination shells in the Heisenberg Hamiltonian. 
In order to ease comparison we used the same exchange 
parameters as in Refs. ITolull 



III. DYNAMICS IN AN EXTERNAL FIELD 
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FIG. 1: Figures (a) and (b) show the distribution of atomic 
moments of the spin dynamics simulations. At finite tem- 
perature the orientations of the atomic spins are distributed 
around a common axis (b). The angles 9 and <j) discussed in 
the text are also defined (a). Figures (c) and (d) show the 
evolution of the spin distribution, as given by the evolution of 
the circular grey disc representing the distribution of magnetic 
moments defined in (b). The system is at finite temperature 
in an external field (c) and in a uniaxial anisotropy (d). 



Different coarse grained levels can be used for describ- 
ing magnetization dynamics. Here we will work with two 
levels: (1) the individual atomic moments, m^, and (2) 
a macro moment, M, representing the sum of the indi- 
vidual atomic moments of the total system. Any atomic 
moment is typically exposed to an interatomic exchange 
field, Bejy, of the order of 1000 T. At finite temperature 
the atomic moments fluctuate around a common direc- 
tion. On average, below the critical temperature, there is 
a finite magnetic moment and the interatomic exchange 
field averaged over all atoms is directed along the average 
moment. The size of the average moment or the macro 
moment depends on the spread of the individual atomic 
moments, which is governed by the temperature. The sit- 
uation is illustrated in the top part of Fig. [TJ where the 
distribution of atomic moments is illustrated for T=0 K 
(Fig. [T^) and at finite temperature (Fig. [TJa) . This de- 
scription of magnetization at finite temperature is the 
starting point for our discussion. 

If the system is exposed to an external field, the aver- 
age moment will precess in this external field. The atomic 
moments precess in a uniform way without distortion of 
their internal distribution. The torque exerted by the 
external field, dnii/dt = — 7111^ x Bgg^ on each atom 
i, results in an equal angular velocity of all the atomic 



moments. This is illustrated in Fig. [3] (top left) where 
the torque (or dmi/dt) is shown as a function of angle 
(9) between the magnetic moment and applied field. In 
Fig. [2] (top right) we also show the resulting angular ve- 
locity (d(j>i/dt) of each atomic spin. The angular velocity 
is constant and seen to be independent of 9, hence the an- 
gular velocity is the same for all spins and it stands clear 
that an external field will not influence the relative orien- 
tation of the atomic spins. The evolution of the distribu- 
tion of the atomic moments at finite temperature during 
relaxation in an external field is schematically shown in 
Fig. QJi. The figure illustrates the fact that an external 
field results in a simple rotation of the magnetization and 
that all individual atomic spins rotate without changing 
the relative direction to all other atomic spins. 



IV. DYNAMICS IN A UNIAXIAL 
ANISOTROPY FIELD 

If there is a uniaxial anisotropy in the system, such as 
magneto-crystalline anisotropy or shape anisotropy, an 
excitation of the macro moment in the anisotropy (by a 
rotation) will in general lead to a precessional motion of 
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the macro moment in the anisotropy field which appears 
similar to the precession in an external field. For the 
case of a uniaxial magneto-crystalline anisotropy which 
we will consider now, there are however important dif- 
ferences in the spin dynamics. We define the anisotropy 
energy for each atomic moment as E — ke 2 z where k is 
the anisotropy constant which determines the strength of 
the anisotropy and e z is the z-component of the direction 
of the atomic moment. The torque and angular velocity 
on any atomic spin are illustrated in Fig. [5] for both an 
easy-axis anisotropy (middle panels) and an easy-plane 
anisotropy (lower panels). The torque is clearly differ- 
ent than in the case of an applied field (top panels), and 
more importantly the angular velocity of each spin is no 
longer independent of 0. Hence, spin dynamics in a uni- 
axial anisotropy field does not lead to a uniform rotation 
of the atomic spins. Instead the internal distribution 
of the atomic moments is distorted during the rotation, 
as illustrated in Fig. [TJi. For the case of the external 
field the distribution of atomic spins remained constant 
with the result that the size of the macro moment re- 
mained constant during the precession. Hence, the pro- 
cess could conveniently be described within a macro mo- 
ment picture. This is not true for the case of the uniax- 
ial anisotropy and the macro moment description breaks 
down. Since the internal distribution of the atomic mo- 
ments is changed during the precession, and the direction 
of any atomic moment in general changes relative to all 
other moments in the system, the size of the macro mo- 
ment changes which leads to a considerably more com- 
plex macro level behavior. In the rotating frame of the 
average moment, the easy-axis anisotropy is seen to coun- 
teract the precession of atomic moments in the effective 
exchange field, while the hard-axis anisotropy is seen to 
enhance the precession in the effective exchange field. As 
the average moment precesses in the uniaxial anisotropy 
the atomic moments will have a tendency to spread re- 
ducing the net moment of the system, as shown in Fig.[5Ji- 
We will refer to this behavior as a spin-wave instability 
(SWI), according to the discussion by Kashuba (Ref. |9|). 

As we will show in our simulations, the SWI results 
in an apparent damping of the uniform motion of the 
macro moment. We define the anisotropy axis as the 
z-axis. What is significant for this damping is the van- 
ishing of the macro moment components perpendicular 
(x, y) to the anisotropy axis and the constant value of 
the parallel macro moment component (z). Hence, dur- 
ing the SWI, the average magnetization of the system 
drops and only the z-component of the average magneti- 
zation remains finite as the x and y-components vanish. 
This gives the relaxation of the macro moment due to the 
SWI a Bloch-Bloembergen form where \M Z \ is constant. 
Thus, due to the SWI there is an alignment of the macro 
moment with the anisotropy axis where the alignment 
occurs maintaining a constant value of the z-component 
of the macro moment. This is illustrated in Fig. [3J This 
shows that a redistribution of angular momentum and en- 
ergy within the magnetic system is taking place. Hence, 
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FIG. 2: The plots illustrate the change of the magnetic mo- 
ment due to an external field, easy-axis anisotropy and an 
easy-plane anisotropy. The graphs on the left hand side give 
the magnitude \dm/dt\ while the graphs on the right hand 
side give the angular velocity of the atomic spins with respect 
to angle 6 between spin and applied field or anisotropy axis. 
Note that in the case of a uniaxial anisotropy field 6 is defined 
as the angle between moment and a fixed crystallographic di- 
rection of the anisotropy field (e.g. 100). H is the strength of 
the external field and K the strength of the anisotropy field. 

there is a relaxation taking place even though the dissi- 
pative damping, a, is set to zero. In reality there is also a 
finite dissipative damping, a, and therefore also a Gilbert 
contribution to the relaxation of the macro moment. In 
some of our simulations, in order to clearly observe the 
Bloch-Bloembergen damping with M z =constant, we set 
a = 0. The fact that the value of the z-component of the 
macro moment is constant during the SWI is expected 
since with zero damping the precessional torque of the 
uniaxial anisotropy is the only source or drain of angular 
momentum within the spin system and this torque lacks 
z-component. 



A. Simulating bec Fe with different strengths of 
uniaxial anisotropy 

In order to study the SWI of bec Fe we choose a 
20x20x20 cell with periodic boundary conditions, en- 
compassing 16000 atomic spins, and three different values 
of the strength of an uniaxial anisotropy: -2 mRy/atom, 
-0.2 mRy/atom and -0.02 mRy/atom, with an easy axis 
directed along the z-axis. Materials with Fe atoms in a 
bec environment and enhanced anisotropy may be found 
experimentally in magnetic multilayers, e.g. with Pt. 
The anisotropy can here be significantly stronger than 
in the bulk case. The magnetic anisotropy of a tetrago- 
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a. Bloch-Bloembergen 
damping 
(|M 2 |=const.) 



b. Gilbert damping 
(|M|=const.) 





FIG. 3: The figures illustrate the Bloch-Bloembergen (a 
damping and the Gilbert (b) damping for a macro moment. 
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FIG. 4: Calculated evolution of the total magnetization of bcc 
Fe as a function of time for different angles between the initial 
magnetization and the applied field and for different values of 
the damping parameter. The temperature was 100 K. 



nal FeCo/Pt(001) superlattice was measured^ to K u — 
2.28 MJm -3 , corresponding to K u « 0.012 mRy/atom. 
The perpendicular magnetic anisotropy of (Co, Fe)/Pt 
multilayers was measured by Sato et alM- to K u — 
0.25 erg/cm -2 , corresponding to K u » 0.027 mRy/atom. 
The strongest magnetic anisotropy found in experiments 
is for SmCdgi&Sa with values of K u = 7.7 MJm- 3 , 
corresponding to K u ks 0.31 mRy/formula unit. As 
will be presented below we see SWI phenomena in our 
simulations for the anisotropy values -2 mRy/atom, - 
0.2 mRy/atom but not for -0.02 mRy/atom. In order 
to simplify, we have not considered non-magnetic (e.g. 
Pt) atoms in the simulations, but only the effect they 
have on the uniaxial anisotropy field. As we will show, 
these systems can display an instability on a time scale of 
picoseconds (shown in Fig. 3]). We now investigate the de- 
pendence of the SWI on thermal fluctuations and damp- 
ing and we investigate the redistribution of the atomic 
moments which takes place. 

In Fig. QJ we show a series of simulations for three dif- 
ferent damping parameters, a, and two different initial 
angles, 9 — 45° and 9 = 90°. In these simulations we 



used a uniaxial energy of -2 mRy/atom. For the macro 
moment there is now a Bloch-Bloembergen like damp- 
ing due to the SWI and a Gilbert like damping due to 
the inclusion of a dissipative damping in the microscopic 
equations of motion. For the case of 9 = 45° we see 
the presence of both these damping terms (see Fig. 0} . 
For a=0.1 the Gilbert term is seen to dominate. After a 
short dip in the magnitude of the magnetization due to 
the SWI the magnitude of the magnetization is seen to 
recover. For cv=0.01 and 0.001 the Bloch-Bloembergen 
damping is seen to dominate, and the size of the mag- 
netic moment reaches a value of M/Mo « 0.6-0.7. For 
9 = 90° the situation is slightly different. At this specific 
angle only the SWI contributes to the relaxation of the 
system. For this reason the behavior in Fig. 3] is fairly in- 
dependent of the magnitude of a, and the magnetization 
evolves with time to a value where M/Mq w 0-0.1. 

The cause of the SWI is an internal redistribution of 
the atomic moments. In Fig. [5] we show a histogram of 
the angles of the atomic moments with respect to the 
average atomic moment. The distribution is shown for 
different points in time for two damping parameters, a = 
0.0 and a = 0.1 (with a uniaxial energy of -2 mRy/atom). 
As a first observation, in contrast to what one might 
expect, the size of a does not change the rate at which the 
directions of the atomic moments are redistributed. This 
is illustrated in both the upper and lower panel of Fig. [5J 
One would, simplemindedly, expect a large damping of 
the atomic moments in the interatomic exchange field to 
reduce the spread of the atomic moments, which would 
counteract the SWI. However, this does not happen. A 
second observation (see upper panel of Fig. [5]) is that the 
distribution of 9 is smeared out during the SWI. This 
is consistent with the fact that the net moment of the 
system is reduced. A third observation (see lower panel 
of Fig. [5]) is that the distribution of <fi is heavily distorted 
during the SWI. At t = the distribution is constant, 
which also is illustrated by the circular disc in Fig.QJl. At 
t = 0.2 ps the distribution is distorted, which is consistent 
with the development of an elliptically shaped disc in 
Fig.[TJl. 

In order to explain the observations in Fig. [5] we show 
in Fig. [6] a histogram of the energy distribution of the 
magnetic moments at different points in time during the 
simulation. The histograms for the energy distribution at 
different points in time fall on top of each other and co- 
incide with the Boltzmann distribution at 300 K, which 
demonstrates that the simulations are done at thermal 
equilibrium, throughout the SWI. This explains the first 
observation of Fig. [5J The effect of the dissipative damp- 
ing in Langevin dynamics is to bring the system to ther- 
mal equilibrium. But since the SWI conserves the ther- 
mal distribution of the system, damping has no net ef- 
fect on the distribution of the directions of the atomic 
moments. The second and third observation from Fig. [5l 
concern the change in angular distribution of the atomic 
moments and explain how the fact that the system re- 
mains in thermal equilibrium can be consistent with a 
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FIG. 5: Distribution of the angles between the average macro 
moment and each atomic moment of the bcc Fe simulation 
cell. In the simulation the initial angle between the average 
magnetization and the anisotropy axis is 9 = 90°. The top 
panel shows the distribution of 9 for the different atomic spins 
and the bottom panel shows the distribution of <j>, defined in 
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reduction in the average magnetization. The angular 
distribution of the atomic moments is heavily distorted 
whereas the energy distribution remains constant. For 
finite damping, the situation changes slightly. We show 
in Fig. [7] how the magnetic energy, which here is the sum 
of exchange and anisotropy energy, evolves in time. For 
the zero damping case in of Fig. [7] the lowering of the 
anisotropy energy is compensated by an increase in the 
exchange energy, leaving the total energy constant. This 
is contrasted by the finite damping cases with a — 0.001 
respective 0.1, were the initial increase in exchange en- 
ergy decays towards its equilibrium value at the given 
temperature. In both cases the time-evolution of the 
atomic moments lowers the total energy. The reason 
for the different behaviors can be explained as follows. 
For the zero damping case, a — 0.000, the sum of the ex- 
change and the anisotropy energy is a constant of motion. 
At the start of the simulations the magnetic moments are 
in thermal equilibrium at T — 300 K and the total ex- 
change energy is constant. When the anisotropy field is 
'turned on' at t = 0, the system is not in an anisotropy 
energy minima as the average magnetization is at an an- 
gle 8 = 90° to the easy axis. The evolving magnetic mo- 
ments lower their anisotropy energy with an amount of 
energy that is in its entity transfered to exchange energy, 





1 1 1 1 1 1 ' 
o o ps 




a -a 0.1 ps 




^^>0.2ps 




&^0.3ps 




— Boltzman distribution 300 K 







10 15 
Energy (mRy) 



FIG. 6: Histogram of the energies of the atomic spins for a 
simulation of bcc Fe with a = 0.0 and 9 = 90°. Although 
there is a large drop of the average moment of the system 
the energy distribution does not change significantly during 
the development of the SWI. Data for different times of the 
simulation are shown. 
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FIG. 7: (color online) The evolution in time of the total en- 
ergy, the exchange energy and the anisotropy energy, for var- 
ious damping parameters. 



since energy can not dissipate in or out of the system. 
With a small but finite damping, a = 0.001, the mag- 
netic excitations can dissipate and lower the exchange 
energy. With a large damping of a = 0.100 the exchange 
energy dissipates to within 5 ps to reach its equilibrium 
value at temperature T — 300 K. 

Thermal fluctuations play an important role for the 
development of the SWI. Naturally there is therefore a 
dependence of the time-scale of the instability on the tem- 
perature. We found however that in the range 10-300 K 
the time-scale is fairly independent on the temperature 
as shown in Fig. [S] (again we used a uniaxial energy of 
-2 mRy/atom for these simulated data). The thermal 
fluctuations also have another role. For systems where 
the macro moment is unable to relax along an anisotropy 
axis (i.e. when — 90°) thermal fluctuations turn out 
as the only mechanism for the system to come out of 
the chaotic SWI state when the anisotropy field is re- 
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FIG. 8: Simulations at different temperatures of bcc Fe with 
a — 0.0 and 6 — 90°. The SWI develops on the same time 
scale for different temperatures. 
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FIG. 9: Starting from bcc Fe in a SWI state the anisotropy 
field is removed (at T=300 K). The system is seen to evolve 
back slowly toward a ferromagnetic state. 

moved. Starting from a chaotic state where the SWI has 
been allowed to bring the system to a zero total moment 
state we suddenly remove the anisotropy field and ob- 
serve the evolution of the system (see Fig. [9]). It is now 
only the complete randomness of the thermal fluctua- 
tions that eventually is able to evolve the system back 
to a ferromagnetic state. The thermal fluctuations will 
eventually bring the spin distribution which has a total 
moment close to zero, to a spin distribution with a to- 
tal moment approaching a finite value. The process is 
however very time consuming, as shown in Fig. [5] and 
only observed for the largest damping parameter in the 
present simulations. 

We now compare simulated results using different 
strengths of the uniaxial anisotropy as well as different 
values of the damping parameter. The interatomic ex- 
change interactions and the size of the simulation cell 
were kept the same as in previous simulations. The 
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FIG. 10: (color online) Calculated evolution of the magneti- 
zation of bcc Fe as a function of time with different values of 
the uniaxial anisotropy and for different values of the damp- 
ing parameter. The magnetization is initially at angle 9 = 90° 
to the anisotropy axis. 



temperature was 300 K. In Fig. [TO] we show the case 
where 0=90, for three values of uniaxial anisotropy, 
—2 mRy/atom, —0.2 mRy/atom and —0.02 mRy/atom. 
Note that the case with an anisotropy of —2 mRy/atom 
was also considered in Fig. [|J although in Fig. [TO] we 
show the dynamical response over a larger time interval, 
5 ps. For the strongest value of the uniaxial anisotropy 
the SWI develops rather easily, whereas for the lowest 
value of the uniaxial anisotropy the SWI does not de- 
velop at all, at least not in the time interval considered. 
The intermediate value of the uniaxial anisotropy results 
in an intermediate situation where the macro moment 
oscillates in time (at least in this time-interval, we will 
return to this situation below). The reason behind the 
different behaviors shown in Fig. [TO] is a competition be- 
tween the strength of the uniaxial anisotropy, which in 
line with the discussion around Fig. [1] tends to spread the 
distribution of all atomic moments, and the importance 
of the other relevant interactions in the system, primarily 
the strength of the interatomic exchange interaction. 

In Fig. [11] we show very similar simulations as in 
Fig. [TOJ with the only difference being that we show re- 
sults for the case when 8 — 45°. Here the intermediate 
and lowest value of the uniaxial anisotropy does not have 
sufficient strength to drive a SWI, whereas the largest 
value of the anisotropy the SWI develops and a Bloch- 
Bloembergen damping occurs. This was also illustrated 
in Fig. 21 but over a shorter time-interval. 

The case when 8 = 45° and with a uniaxial anisotropy 
of -0.2 mRy/atom is, as Fig. [TO] suggests, a particularly 
interesting case, since here the anisotropy and exchange 
interactions seems to be tuned into a situation where 
both are very influential for the evolution of the macro 
spin. In fact, Fig. [TOI suggests that in this case the mag- 
netization oscillates between a Gilbert like damping and 
Bloch-Bloembergen like damping. For this reason we 
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FIG. 11: (color online) Same as Fig. [I0]but with the magne- 
tization initially at angle 9 — 45° to the anisotropy axis. 



have extended the simulations over a larger time interval 
(50 ps), and the results are shown in Fig. [T3J It is to 
be noted from this figure that for this borderline case, 
the evolution of the macro spin depends not only on the 
competition between interatomic exchange and uniaxial 
anisotropy, but also on the value of the damping param- 
eter. For large values of the damping a regular Gilbert 
damping behaviors is found. For small and intermediate 
values of the damping the macro spin is found to oscil- 
late in time, but otherwise following a dynamic response 
which resembles Bloch-Bloembergen damping. Hence the 
data in Fig.[T2lshow that by careful tuning of the relative 
importance of the uniaxial anisotropy, exchange interac- 
tion and damping, one may obtain a behavior which is 
more complex than that given by pure Gilbert or Bloch- 
Bloembergen damping. 

The finite size of the simulation cell restricts the pos- 
sible spin wave excitations. The simulations described so 
far were all for L = 20 corresponding to 16000 magnetic 
moments. With a smaller cell only the modes with short 
wave lengths can occur. This means that the weaker 
uniaxial anisotropy cannot drive an SWI unless the sim- 
ulation cell is large enough. The trend for if a SWI can 
occur or not for the different simulation cells, with cell 
size L = 10, 15, 20 and 25, are presented in Tables HI fTTl 
for the 90 degree and 45 degree case, respectively. In the 
table for the 90 degree case we have defined a strong SWI 
as the case where the magnetization drops below 0.2M , a 
medium SWI as when it drops below 0.6M , a weak SWI 
as when the magnetization drops with 0.05 — 0.20Mo and 
no SWI when the magnetization drops less than 0.05Mo. 
In the table for the 45 degree case the same notation is 
used apart from that we here redefine strong SWI as when 
the magnetization drops below 0.65Mo (which here corre- 
sponds to a Bloch-Bloembergen damping). The results of 
Tables. Hill] correspond well to the results of Ref. [13 (see 
e.g. Eqn. 27). For the anisotropy values— 2 mRy/atom, 
—0.2 mRy/atom and —0.02 mRy/atom and with the ex- 
change energy summed up over all coordination shells to 



TABLE I: Simulations for varying cell size, 90 degree case. 
The entries describe the possible occurrence of a SWI during 
the simulation time t = 5 ps. 



K u (mRy/atom) 


a 


L=10 


L=15 


L=20 


L=25 


-2.0 


0-0.1 


strong 


strong 


strong 


strong 


-0.2 


0-0.1 


no 


weak 


medium 


medium 


-0.02 


0-0.1 


no 


no 


no 


no 



TABLE II: Simulations for varying cell size, 45 degree case. 
The entries describe the possible occurrence of a SWI during 
the simulation time t = 5 ps. 



K u (mRy/atom) 


a 


L=10 


L=15 


L=20 


L=25 


-2.0 


0-0.01 


strong 


strong 


strong 


strong 


-2.0 


0.1 


weak 


weak 


weak 


weak 


-0.2 


0-0.01 


no 


no 


no 


weak 


-0.2 


0.1 


no 


no 


no 


no 


-0.02 


0-0.1 


no 


no 


no 


no 



w 10 mRy/atom we get N max = 5, 16, 51 where N„ 
the largest cell size that suppress SWI effects. 



V. DISCUSSION AND SUMMARY 

In this paper we have investigated the conditions when 
a spin-wave instability (SWI) may occur. In order for 
this to happen, a number of requirements must be met. 
First, there must be an initial perturbation to the system, 
e.g. thermal fluctuations, such that the atomic moments 
start to deviate from the direction of the macro moment. 
Secondly, there must be a magnetic anisotropy in the sys- 
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a=0.1000 
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20 30 
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40 50 



FIG. 12: (color online) Same as the middle panel of Fig. \W\ 
but showing the evolution of the magnetization up to 50 ps. 
For a = . . . 0.01 the magnetization oscillates in the interval 
0.55 — 0.70Mo. For a — 0.1 the magnetization recovers after 
~ 20 ps to the value M = O.88M0 which is the thermally 
equilibrated value at temperature T = 300 K. 
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tern. The presence of a SWI is found to a large degree 
be determined by a competition between the magnetic 
anisotropy and the strength of the exchange interaction. 
In some special cases, where these two contributions are 
very delicately balanced, the value of the damping pa- 
rameter can finally determine whether or not a SWI oc- 
curs. We have also found that the size of the simulation 
cell is influential for if a SWI occurs, a conclusion which 
is in agreement with the results of ReffiS 

Another conclusion we reach from our simulations is 
that due to thermal fluctuations the simple model of 
a macro moment precessing in a uniaxial anisotropy is 
found to be inaccurate. The uniaxial anisotropy leads 
to a non-uniform rotation of the composing atomic mo- 
ments. On a short time scale the effect is small. On 
a longer time scale or for larger anisotropies there are 
severe consequences. An instability appears which ef- 
fectively leads to a Bloch-Bloembergcn damping of the 
magnetization. 

Our simulations point to a technical avenue for de- 
signing media for data-storage and magnetic memories, 
where e.g. the grain size of the storage media would be a 
materials property which one could compare to the vari- 



ous sizes of our simulation cell. Media with a small grain 
size could possibly then exhibit a weaker tendency for 
a SWI to be observed. If experimental evidence for the 
spin wave instability could be demonstrated, it would 
imply that there is an increased importance to a fine 
grain description of the magnetization dynamics in sim- 
ulations and it would show that macro moment mod- 
els lose accuracy when anisotropies are involved in the 
dynamics. Further experimental studies addressing this 
issue are highly desired. 
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